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Objectives 

The  primary  objectives  of  this  work  are  the  construction  of  a  rigorous  mathematical 
framework  and  the  corresponding  computational  science  tools  that  can  be  used  to 
address  problems  of  parameter  identification,  real-time  tracking  and  estimation  for 
spatially  dependent  systems.  This  includes  determining  optimized  sensor /actuator 
locations  for  complex  hybrid  spatial  systems  to  enhance  tracking,  estimation,  infor¬ 
mation  and  effectiveness  while  limiting  energy  consumption.  Reduced-order  modeling 
techniques  are  implemented  as  an  efficient  way  to  compute  the  functional  gains.  We 
illustrate  how  a  small  number  of  strategically  placed  sensor /actuator  is  sufficient  to 
stabilize  the  flow  while  inappropriate  placement  of  these  sensors  could  destabilize 
the  flow.  Additionally  we  consider  information  delays  present  in  the  sensor /actuator 
network.  The  models  are  complex  multi-scale  systems  of  coupled  partial  and  delay 
differential  equations.  We  show  that  under  suitable  conditions,  the  coupled  delay 
PDE  systems  are  well  posed  and  we  use  this  corresponding  abstract  formulation  to 
construct  efficient  numerical  methods  for  control  design. 

Overview 


Sensor  Placement 

Feedback  control  and  optimal  state  estimator  problems  for  PDE  systems  lead  to 
operator  equations  where  the  feedback  and  observer  gain  operators  often  have  integral 
representations.  These  kernel  functions  are  called  functional  gains  and  often  have 
localized  support  that  provides  information  about  the  controller  and  estimator.  If 
these  functional  gains  exist  and  can  be  computed,  one  can  use  this  information  to 
place  sensors  in  those  regions  in  space  that  are  most  important  to  the  controller  and 
estimator. 

The  functional  gain  operators  are  computed  by  solving  appropriate  approximations  to 
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the  infinite  dimensional  Riccati  equations  arising  from  infinite  dimensional  estimation 
and  control  problems.  Due  to  the  large  scale  of  the  approximate  problem,  special 
techniques  are  required  to  solve  these  equations.  We  propose  and  demonstrate  a 
reduced-order  modeling  technique  to  compute  the  functional  gains  and  analyze  these 
gains  to  obtain  insight  into  the  physical  behavior  of  the  feedback  control  mechanism. 

For  the  discussion  consider  a  flow  past  a  circular  cylinder  at  Re=200  and  we  use 
suction  on  the  cylinder  surface  to  suppress  vortex  shedding.  The  incompressible 
Navier-Stokes  equations  are  solved  using  parallel  algorithms  and  snapshots  of  actu¬ 
ated  and  unactuated  flow  data  are  collected  for  building  a  reduced-order  model.  In 
this  case,  the  reducedO-order  model  is  computed  using  proper  orthogonal  decomposi¬ 
tion  (POD). 

The  velocity  field  is  written  as  the  sum  of  the  mean  flow  (v)  and  the  velocity  fluctu¬ 
ations  (v7).  The  fluctuations  are  expanded  in  terms  of  the  POD  eigenfunctions  (<3>j) 
computed  from  the  collected  snapshots  and  the  control  function  method: 

M  m 

v(x,t)  «  v(x)  +  ^z[(t)$i(x)  +  J^ui(t)ri(x).  (1) 

i= 1  i= 1 

The  number  of  POD  modes  is  denoted  by  M  and  m  is  the  total  number  of  control 
modes.  Each  function  rj(x)  is  a  suitable  divergence- free  control  function  that  sat¬ 
isfy  the  inhomogeneous  boundary  condition  due  to  fluidic  actuators,  and  u(t)  is  the 
variable  control  input. 

We  project  the  Navier  Stokes  equations  onto  the  POD  basis  functions  (or  modes)  to 
develop  a  reduced-order  model.  The  reduced-order  model  for  rn  —  1  has  the  form 


=  A(  +  A 
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The  system  is  then  linearized  about  its  mean  component  to  reduce  the  model  to  a 
standard  linear  time  invariant  (LTI)  control  problem.  If  the  control  is  a  linear  function 
of  the  states,  then  the  closed  loop  system  then  takes  the  form 
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As  the  purpose  of  the  control  law  is  to  stabilize  the  system,  it  is  desired  to  have 
the  poles  of  the  system  in  the  left-half  of  the  complex  plane.  In  other  words,  the 
eigenvalues  of  Ac  must  have  negative  real  parts. 

The  feedback  control  gain,  K,  represents  an  approximation  to  the  feedback  operator 
associated  with  the  PDE  control  system.  The  Ricsz  Representation  Theorem  implies 
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that  there  exists  a  divergent  free  vector  held  h(a:)  such  that 


M 


uc(t)  =  -K,z r(t)  -  Kuu(t )  =  -  h(x)  •  ^2  zl(t)$i(x)dQ,  -  Kuu(t).  (4) 


i=  1 


The  efficiency  of  this  approach  was  verihed  in  the  following  way:  For  a  specific  con¬ 
troller  where  the  eigenvalue  is  moved  to  —0,05  ±  1.16*,  12  POD  inodes  were  used 
compute  the  functional  gain,  see  Figure  1  .  From  this  information,  we  simulate  two 


Figure  1:  Functional  gain  approximated  at  four  locations;  good  (square)  and  bad  (circle) 
locations. 

scenarios.  In  each  of  the  scenarios  we  use  four  locations  to  approximate  the  feedback 
gain.  These  choices  are  motivated  by  the  magnitude  of  the  functional  gain  in  the 
wake.  Four  locations  are  chosen  where  the  functional  gain  is  high  (good  choices)  and 
four  where  the  gain  has  small  magnitude  (bad  choice).  These  locations  are  identified 
in  Figure  1  with  good  choice  and  a  bad  choice  scenarios  represented  by  squares  and 
circles,  respectively.  We  approximate  the  feedback  control  input  using  quadrature  as 
follows: 


where  dA  —  dx  ■  dy  and  dx  —  dy  —  0.2.  The  approximation  of  the  feedback  gains 
provided  by  the  good  data  set  stabilizes  the  system.  Although  the  response  is  slow 
compared  to  the  full-order  feedback  control,  it  is  important  to  note  that  the  gains 
were  approximated  only  from  four  locations  rather  than  using  the  complete  domain. 
On  the  other  hand,  if  the  sensors  are  placed  at  bad  locations,  the  feedback  gain 
approximation  is  poor  and  fails  to  stabilize  the  closed-loop  system.  Rather,  for  long¬ 
time  integration,  it  results  in  instability  of  the  system.  Thus,  it  is  important  to  choose 
the  correct  locations  of  the  sensors  to  measure  or  estimate  the  flow  held. 

By  computing  the  functional  gains  and  using  the  spatial  information  contained  in 
these  gains  we  could  determine  where  to  place  sensors  in  the  wake  region  to  ensure 
stability  and  performance.  The  control  is  computed  efficiently  by  using  good  quadra¬ 
ture  points  for  the  integral  operator. 
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Figure  2:  Closed- loop  simulation  of  u  with  Figure  3:  Closed- loop  simulation  of  u  with 
gains  approximated  at  good  locations.  gains  approximated  at  bad  locations. 

Observe  that  this  analysis  also  provides  valuable  information  concerning  the  construc¬ 
tion  of  state  estimators  (filters).  In  particular,  if  one  can  not  physically  place  a  sensor 
in  the  desired  location  in  the  wake,  then  one  would  design  a  filter  to  estimate  the  flow 
only  in  that  region.  Thus,  this  approach  can  be  used  for  designing  reduced-order 
controllers  without  developing  reduced-order  models  first. 

Delayed  actuator  dynamics 

This  work  investigates  the  inclusion  of  actuator  dynamics  with  delays  on  a  control 
problem.  In  many  applications  such  as  the  control  of  energy  efficient  buildings  where 
actuation  is  provided  by  a  HVAC  system,  the  inclusion  of  actuator  dynamics  results 
in  a  more  realistic  representation  of  the  system.  Including  the  actuator  dynamics 
impacts  the  support  of  the  functional  gain  k0(-)  which  in  turn  impacts  sensor  loca¬ 
tion  problems.  Parabolic  partial  differential  equations  with  no  actuator  dynamics  are 
typically  formulated  as  an  abstract  control  problem.  Coupling  the  actuator  dynam¬ 
ics  defined  by  a  delay  equation  with  the  abstract  system  at  this  stage  results  in  a 
composite  system  that  can  be  difficult  to  analyze.  Computational  schemes  based  on 
finite  element  methods  often  require  high  order  “test  functions” . 

We  develop  an  abstract  state  space  model  that  avoids  many  of  these  complexities  by 
working  directly  with  the  PDE  system.  Standard  finite  element  methods  can  be  used 
to  construct  approximations  of  the  composite  control  system  and  in  many  cases  these 
approximations  can  be  implemented  with  existing  software  packages. 

Consider  the  convection  diffusion  control  system 

0  =  /iV20(t,  0  +  k(£)  ■  V0(t,  0,  (5) 

with  Dirichlet  boundary  inputs  9(t,  £)|an  =  g(0vo(t)i  where  g(-)  G  L2(dfl)  is  a  fixed 
function  defined  on  the  boundary.  In  this  setting  g(-)  can  have  compact  support,  in 
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which  case,  the  input  to  the  PDE  (5)  is  local.  Following  the  development  in  [2],  let 
b(-)  be  the  weak  solution  of  the  boundary  value  problem 

/iV26(0  +  «(0  •  V6(£)  =  0,  6(0 |an  =  <?(0- 

We  are  concerned  in  the  case  where  the  input  vg(t)  is  given  by  the  output  of  a  finite 
dimensional  delay  differential  equation  describing  the  actuator  dynamics.  Assume 
that  vg(t)  =  Hw(f) 

where 

w (t)  =  A0w (t)  +  A]w(t  -  r)  +  Bau(f).  (6) 

ffere,  A0  and  A]  are  nxn  matrices,  Ba  is  an  n  x  m  matrix  and  r  >  0  is  a  fixed  time 
delay,  and  u(-)  :  [0,+oo)  — y  Mm  is  the  command  to  the  actuator.  Assuming  one  can 
differentiate  vg(-)  it  follows  that  vg(t)  =  HA0w(f)  +  HAiw(f  —  r)  +  HB„u(f).  Let 
9(t,  0  =  0  (t,  0  —  b(£)vg(t)  then  the  composite  system, 

^0  ( t ,  0  =  [hV2  +  «(0  ■  V]  9  ( t ,  0  +  f0(z)w (t)  +  fi(x)w(f  -  r)  +  be(x)u(t),  (7) 

w (t)  =  A0w (t)  +  Apw (t  -  r)  +  Bau(f)  (8) 

can  be  written  as  a  distributed  parameter  delay  system  on  the  state  space  X  = 
L2(fl)  x  Rn 

x(t)  =  Aox(t)  +  A\x{t  —  r)  +  B0u(t).  (9) 

The  following  result  follows.  See  [9]  for  the  details  on  the  proof  and  notation  which 
are  omitted  due  to  lack  of  space. 

Theorem  1  The  linear  operator  Aq  generates  a  Co-semigroup  T0(t )  on  X  =  L2(Tl)  x 
Mn.  Moreover,  if  Ag  and  A0  are  stable  and  crp(Ae)ncrp(A0)  =  0,  then  the  composite 
semigroup  is  exponentially  stable  and  the  controlled  system  (9)  is  stabilizable. 

Let  z(t)  =  [  9(t ,  •)  w(f)  wt(-)  ]  =  [  9(t ,  ■)  <3>(t)  ]  then  one  can  also  show,  see 

[10]  for  the  detail,  the  composite  distributed  parameter  control  system  has  the  form 

z(t)  =  Acz{t )  +  Bcu{t)  EZ  =  L2(fl)  x  pET  x  L2(-r,  0;  Mn)],  (10) 

with  bounded  (compact)  input  operator.  The  following  result  follows,  see  [10]  for  the 
detail. 

Theorem  2  The  linear  operator  Ac  generates  a  Co-semigroup  Tc(t )  on  Z  =  L2{Vt)  x 
Mn  x  L2(— r,  0;Mn).  Moreover,  if  Ag  and  Ag  are  stable  and  crp(Ag)r\ap(Ad)  =  0, 
then  the  composite  semigroup  is  exponentially  stable  and  the  controlled  system  (10) 
is  stabilizable. 
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Consider  the  Linear  Quadratic  Regulator  (LQR)  problem  for  the  composite  system 
(10)  with  cost 


/*+oo  /»+oo 

Jc=  {{Qz(t),z(t))z}+  /  {(Rau(t),  u(f))Rm}  dt.  (11) 

J  0  Jo 

The  details  of  this  transformation  are  described  in  [9]  and  will  not  be  presented  here 
in  order  to  save  space.  Implementing  the  Riesz  representation  theorem  results  in  the 

feedback  controller  of  the  form  uop*(f)  =  —  J  k e{^)0(t,f)df  —  kaw(f)  —  kd(s)w(f  + 

n 

s)ds,  where  k“  =  ka  —  kd(£)b(£)Hdf.  The  functions  fcf(-)  and  kfj(-)  are  called  the 

functional  gains  and  ka  and  kJ(£)  are  mxn  matrix  functions.  Consequently,  once  one 
computes  ke(£)  and  ka,  computing  the  gain  matrix  ka  requires  only  a  quadrature. 

Numerical  Example 

We  consider  the  ID  control  problem  for  the  parabolic  PDE 

^  e  (0, 1),  and  0(t,  0)  =  v0(t),  G(t,  1)  =  0.  (12) 


Use  the  following  parameters:  ji 
function  b(-)  is  given  by  b(x)  = 

Dip(-)  =  J1  d(x)cp(x)dx  and  d(x) 


(e 


1/120,  k  =  .1,  R 
kx  -  e~k) / (1  -  e~k) 
f  5,  0.4  <  x  <  0.6 
1  0,  elsewhere 


=  0.1,  Ra  =  0.1,  and  the 
.  Set  Q  =  D*  x  D  where 

.  The  actuator  is  modeled 


by  the  simple  scalar  delay  equation  w(t)  =  —1.0 w(t)  +  0.95 w(t  —  r)  +u(t )  and  vo(t)  = 
Htu(t)  =  5 w(t).  A  standard  piecewise  linear  finite  element  scheme  (with  upwinding) 
is  used  to  approximate  the  convection-diffusion  equation  (see  [6],  [' ]  and  [8])  and 
the  “AVE”  scheme  is  used  to  approximate  the  delay  system  (see  [1]  for  details). 
Subdivide  the  interval  (0,1)  into  N  equal  subintervals  and  the  interval  (— r,  0)  into 
M  subintervals.  Due  to  space  limitations,  we  will  not  provide  a  rigorous  proof  of 
convergence  here.  However,  the  numerical  results  presented  below  illustrate  this 
convergence. 


Figure  4  illustrates  the  convergence  of  the  functional  gains  k^M(-)  and  k^M(-).  We 
note  that  the  system  (10)  is  stable  since  maxQi  =  max(real(Ac ))  =  —0.0255.  Also, 
the  closed  loop  operator  satisfies  maxci  =  max{real([Ac  —  JCBc]))  =  —0.1592. 


The  inclusion  of  actuators  offers  a  more  “realistic”  model  of  the  dynamics  of  the 
coupled  system  and  allows  the  integration  of  system  components.  Including  the  actu¬ 
ator  dynamics  impacts  the  support  of  the  functional  gain  k0(-)  which  impacts  sensor 
location  problems  (see  [3],  [4],  [5]),  and  [12].)  Similar  results  were  observed  for  the 
delayed  actuator  case,  see  [9]  and  [10]  for  more  detail. 
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(a) 


(b) 


Figure  4:  Plots  of  (a)  ke n,m(')  and  (b)  kdN,M{--) 


Infinite  Dimensional  Delay  Differential  Equations  in  Control  and 

Sensitivity  Analysis 


Here  we  also  consider  a  class  of  functional  partial  differential  equations  (FPDEs) 
which  arise  naturally  in  problems  of  control  of  systems  governed  by  partial  differen¬ 
tial  equations  where  delayed  actuator  dynamics  are  included  and  in  the  sensitivity 
analysis  of  such  systems  when  one  is  concerned  with  sensitivities  with  respect  to  de¬ 
lays.  Numerous  examples  motivating  these  models  are  presented  in  [  ].  These  exam¬ 

ples  include  Boundary  Control  of  PDE  Systems  with  Actuator  Dynamics,  Boundary 
Control  with  Delayed  Actuator  Dynamics,  Distributed  Control  with  Neutral  Delayed 
Actuator  Dynamics,  Boundary  Control  with  Neutral  Delayed  Actuator  Dynamics, 
and  Sensitivity  Systems.  Well-posedness  for  a  class  of  FPDE  systems  in  product 
spaces  is  established  and  these  formulations  are  used  as  a  framework  to  develop  ef¬ 
ficient  numerical  approximations  for  control  and  simulation  of  the  PDE  problems. 
The  theoretical  results  extend  existing  well-posedness  results  to  problems  where  the 
standard  range  condition  does  not  apply  and  we  present  a  conjecture  about  a  more 
general  theorem. 

We  consider  the  control  systems  defined  by  a  class  of  infinite  dimensional  neutral 
functional  differential  equations  (NFDEs)  of  the  form 


z(t)  —  D z(t  —  r)  =  A oz(t)  +  Aiz(t  —  r )  +  B u(t)  €  Z  =  0  xR" 
where  0  is  a  Hilbert  space  and  D,  Ao  and  Ai  have  the  specific  structure 

Be  ' 

B  J  ’ 

respectively. 


0  D 

0  E 


— 


Ag  F0 
o  Ao 


u  — 


0  F1 

0  Ai 


and  B  = 


(13) 


(14) 
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When  the  operator  D  is  not  zero,  well-posedness  for  a  generic  system  of  the  form  (13) 
often  requires  assuming  the  “range  condition” 

Range( D)  C  Domain^ A0)  and  A0O  is  bounded.  (15) 

ffowever,  for  certain  applications  to  control  and  sensitivity  analysis  of  PDE  systems 
with  delays,  the  range  condition  (15)  fails.  In  particular,  boundary  control  and  sen¬ 
sitivity  analysis  of  PDE  systems  where  delayed  actuator  dynamics  are  included  can 
lead  to  systems  of  the  form  (13)  where  (15)  is  not  satisfied.  We  formulate  the  system 
(13)  defined  by  (14)  as  an  abstract  system  on  a  “reduced”  product  space  and  take 
advantage  of  the  problem  structure  to  obtain  well-posedness  without  the  range  con¬ 
dition  (15)  and  use  this  framework  to  develop  efficient  approximations.  We  assume 
the  following  hypothesis  holds. 


Hypothesis  H 

(HI)  The  linear  operator  Ae  :  D(Ag)  C  0  — y  0  generates  a  Co-semigroup 
semigroup  on  0. 

(H2)  The  linear  operators  D  :  Mn  — >  0,  Be  :  — >  0  and  Ft  :  Mn  — )■  0, 

i  —  0, 1,  are  bounded. 

(H3)  The  finite  dimensional  operators  E  :  — >■  Mn,  A*  :  Mn  — >  Mn, 

i  —  0, 1,  and  B  :  Mm  — x  M.n  are  linear  (hence  bounded). 


The  state  z(-)  G  Z  =  0  x Mn  has  the  form  z(t )  =  [  9(t,  •)  w(t)  ]T  Z  =  <3  xRn  and 
because  of  the  structure  (14)  only  the  finite  dimensional  state  w(-)  contains  delays. 
With  this  in  mind  we  define  the  (reduced)  product  space 


XR  =  Z  x  L2((— r,  O);P)  =  0xrx  £2((-r,  0);  W1) 

and  note  that  XR  may  be  viewed  as  a  subspace  of  X  =  Z  x  L2((— r,  0);  Z)  by  the  nat¬ 
ural  embedding.  Consider  the  abstract  control  system  on  XR  =  Z  x  L2((— r,  0);  Mn)  = 
0xPx  L2((— r,  0);  Mn)  defined  by 

x(t )  =  Ax(t)  +  Bu(t)  G  XR,  (16) 


where  the  domain  of  A  is  given  by 
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D{A)  = 


n 

¥>(•) 


G  D{Ao)  x  Mn  x  H\[-r,  0];  Mn)  :  y>(0)  -  E( p(-r)  =  rj  V  ,  (17) 


and  for  x  =  [  9  rj  <p(-)  ]T  G  D(A) 


e 

Ae6  +  F0<p(0)  +  Fitp(-r) 

V 

= 

A0^(0)  +  Ai  <p(-r) 

.  ^(-) . 

i 

■«i-a 

Ax  =  A 


(18) 


The  following  result  applies  to  the  special  case  where  (H 4)  holds  instead  of  the  range 
condition  on  D. 

Theorem  3  Assume  (HI)  -  (H 3)  hold.  If  either 

(HD)  the  operator  D  satisfies  the  the  reduced  range  condition 

Range(D)  C  Domain(Ag)  and  AqD  is  bounded  or 

(HA)  the  bounded  linear  operators  Fq  and  F\  are  given  by 

F0  =  FA0  and  F\  =  FAi  (19) 

then  the  operator  A  defined  by  (17)  -  (18)  generates  a  Co-semigroup  X(t)  on  XB  = 
0  x  Rn  x  L2((-r,0);Mn). 

The  detail  and  proof  of  this  result  is  presented  in  [11], 

Numrical  example 

We  use  the  formulation  (16)  to  construct  approximating  systems  for  simulation  and 
control.  For  example,  we  apply  this  approach  to  the  parabolic  boundary  control  prob¬ 
lem  with  actuator  dynamics  described  by  a  (retarded)  delayed  differential  equation. 
In  one  dimension  this  problem  comes  from  the  boundary  control  problem  defined  by 
the  parabolic  PDE 

r\  r\  2  r\ 

gjT(t,0^^T(t,0  +  ^T(t,0,  C  e  (0,1),  (20) 

with  boundary  control  at  the  left  end  and  Dirichlet  condition  at  f  =  1  so  that 

T(t,0)  =  v(t),  T(t,  1)  =  0.  (21) 

Assume  that  the  actuator  dynamics  are  given  by  the  finite  dimensional  retarded  delay 
differential  equation 


w(t)  =  A0w(t)  +  A  iw(t  —  r)  +  B  u(t),  (22) 

where  r  >  0.  Also,  for  simplicity  we  assume 

v(t)  =  Hw(i  —  r).  (23) 

Consider  an  LQR  problem  defined  by  the  boundary  control  problem  (20)  -  (21)  with 
cost  function 


J 


{(QT(t,  -),T(t,  -))L2  +  (Qaw(t),w(t))}  dt, 


9 


where  Qa  =  [H]TAH  >  0. 

For  the  numerical  computations  we  use  the  parameters  /i  =  1/240, 
k  =  .01,  R  =  .1  and  set  Q  =  D*D  where  D  :  L2( 0, 1)  — y  M1  is  defined  by 


D<p(-)  =  /  d(x)ip(x)dx 


and  d(-)  e  L2( 0, 1)  is  given  by 


f  5,  0.4  <  x  <  0.6 
1  0,  elsewhere 


Also,  observe  that  since  k  —  -  —  24  >  0,  the  function  b(-)  is  given  by  b(x)  = 
(e ~kx  —  e~k)/(l  —  e~k).  The  actuator  is  modeled  by  the  simple  scalar  delay  equation 


w(t)  =  —1.25  w(t)  +  1.0  w(t  —  r)  +  u(t). 


The  actuator  dynamics  are  assumed  to  be  given  by  (22)  and 


Veit)  =  0.5  w(t  —  1) 


so  that  H  =  0.5.  The  weight  on  the  input  v(-)  is  set  to  R  =  0.1  and  the  actuator 
control  weight  is  set  to  Ra  =  0.1.  A  standard  piecewise  linear  finite  element  scheme 
(with  upwinding)  is  used  to  approximate  the  convection-diffusion  equation. 

Again,  it  is  important  to  note  that  only  w{-)  is  delayed  which  implies  that  the  feedback 
control  law  does  not  require  the  past  history  of  9(t,  •)  so  that 


u opt(t) 


j  ke(£)9(t,  £)dx  —  kaw(t) 
n 


kd(s)w(t  +  s)ds 


and  the  goal  is  to  compute  the  functions  ke[-)  and  kd[-)  and  a  constant  gain  ka. 
A  second  implication  of  the  fact  that  only  w(-)  is  delayed  is  that  this  allows  us 
to  use  a  reduced  finite  volume  method  (the  so  called  “reduced  AVE”  scheme)  for 
approximating  the  composite  delay  system. 


We  subdivided  the  interval  (0, 1)  into  N  equal  subintervals  for  N  =  16,  32,  64, 128,  256, 
512, 1024  and  the  interval  (— r,  0)  into  M  equal  subintervals  for  M  =  32,  64, 128,  256, 
512, 1024,  2048.  A  standard  finite  element  approximation  of  the  convection  diffusion 
equation  on  (0, 1)  is  then  coupled  to  a  piecewise  constant  approximation  of  the  past 
history  terms  on  (— r,  0).  Note  that  we  used  M  =  2N  since  the  finite  element  scheme 
has  a  convergence  rate  twice  the  AVE  scheme.  We  solve  the  algebraic  Riccati  equation 
using  Matlab  and  compute  the  functional  gains  and  gain  parameters  that  define  the 
optimal  LQR  feedback  controller. 
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N 

M 

ka 

16 

32 

0.9897 

32 

64 

0.9812 

64 

128 

1.0240 

128 

256 

1.0726 

256 

512 

1.0994 

512 

1024 

1.1115 

1024 

2048 

1.1118 

Tabic  1:  Values  of  k%M 


Figure  5  illustrates  the  convergence  of  the  functional  gains  k9  and  in  Figure  6 

we  see  the  same  for  the  functional  gains  kdN,M(-)-  Note  that  the  functional  gains 
ke n,m{')  f°r  the  convection  diffusion  equation  converge  at  N  =  64  so  we  only  show 
plots  of  k9 n,m(')  f°r  N  =  16,32,64  and  M  =  32,64,128.  As  shown  in  Figure  6, 
convergence  of  the  delay  functional  gains  is  much  slower.  Finally,  Table  1  contains 
values  of  ka  as  N  and  M  increase. 


Figure  5:  Plots  of  ke n,m(') 

The  reduced  AVE  scheme  for  the  delayed  terms  was  selected  for  two  reasons.  First, 
this  scheme  meets  all  the  conditions  necessary  to  establish  the  type  of  convergence 
required  for  gain  convergence.  The  second  reason  is  that  the  reduced  scheme  greatly 
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Figure  6:  Plots  of  kdN:M(-) 

“reduces”  the  overall  size  of  the  approximation  scheme.  For  example,  if  we  had 
used  the  full  AVE  scheme  where  N  —  1,  024  and  M  =  2,  048,  then  the  size  of  the 
approximating  system  would  have  been  N(M +1)  =  2, 098, 176.  However,  the  reduced 
system  has  is  N  +  M  =  3,  072  so  there  is  a  huge  reduction  in  model  size.  Finally,  to 
obtain  reasonable  convergence  for  this  specific  example,  we  needed  to  use  M  =  2,  048 
approximations  for  the  delay  system. 

We  establish  the  well-posedness  for  a  class  of  these  systems  in  product  spaces  and 
use  this  formulation  as  a  framework  to  develop  efficient  numerical  approximations 
for  control  and  simulation  of  the  PDE  control  problems. 
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